Out-of-plane spin polarization from in-plane electric and magnetic fields 
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We show that the joint effect of spin-orbit and magnetic fields leads to a spin polarization per- 
pendicular to the plane of a two-dimensional electron system with Rashba spin-orbit coupling and 
in-plane parallel dc magnetic and electric fields, for angle-dependent impurity scattering or non- 
parabolic energy spectrum, while only in-plane polarization persists for simplified models. We 
derive Bloch equations, describing the main features of recent experiments, including the magnetic 
field dependence of static and dynamic responses. 



Generating spin populations at a nanometer scale is 
one of the central goals of spintronics Using spin- 
orbit interaction promises electrical control, allowing to 
integrate spin generation and manipulation into the tra- 
ditional architecture of electronic devices. Bulk spin po- 
larization, driven by electron drift in an electric field, 
was predicted long ago for noncentrosymmetric three- 
(3D) and two-dimensional (2D) systems 0, S, S H B • 
In 2D, the polarization is in-plane, typically along the 
effective spin-orbit field bdr — (bso(k)} ^ 0, obtained 
by averaging spin-orbit coupling over the distribution of 
electron momenta fi-k ||8|. In-plane polarization compo- 
nents were observed recently in p-GaAs hetero junctions 
Pi , quantum wells T^l , and strained n-InGaAs films . 
Out-of-plane spin polarization can be generated by the 
spin-Hall effect, but only near sample edges Be- 
low, we propose a mechanism for out-of-plane spin po- 
larization generated in the bulk by applying an in-plane 
magnetic field B. This perpendicular polarization allows 
efficient optical access, e.g., via Kerr rotation. We find 
that the use of such an average field bdr is not always 
valid. Naively, one might consider the system as being 
subject to a total in-plane field (b), given by the sum 
of B and bdr, see Fig. 1(a). In steady state, one then 
expects electrons to be polarized along this total field: in 
particular, no polarization perpendicular to the bdr, B 
plane. Algebraic addition of these fields worked well in 
describing Hanle precession of optically oriented 2D elec- 
trons in GaAs jl3||. However, Kato et al. ^] reported a 
spin polarization that is incompatible with such a naive 
picture and emphasized the need of identifying its micro- 
scopic mechanisms. 

In this article, we develop a theory describing the inter- 
play between spin-orbit interaction and external electric 
and magnetic fields in the presence of impurity scatter- 
ing, and demonstrate that the concept of average spin- 
orbit field is subject to severe restrictions. The naive 
expectation turns out to be correct only in the special 
case of parabolic bands and isotropic impurity scatter- 
ing. However, as we show below, for anisotropic scatter- 
ing (e.g., small angle scattering), such correlations result 
in a more complex structure of the distribution function 
and an out-of-plane spin polarization. Concretely, in- 
terplay of bdr o.'Tid B leads to a generation term in the 
Bloch equation proportional to bdr x B whose magnitude 
is controlled by anisotropy of potential scattering and 
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Figure 1: (color) (a) Field geometry, assuming g/iB > 0, a < 
0. Out-of-plane spin polarization is electrically generated with 
rate F^ (blue arrow) due to the interplay between spin-orbit 
interaction, external electric field E and magnetic field B, and 
anisotropic impurity scattering. The polarization precesses 
(blue arc) in p'/^bB -|- (bso)- (b) Dynamics of out-of-plane 
component of spin polarization generated by a short electrical 
pulse of length tp < 1 ns, for 70 = 0.03, g* = 0.65 and = 
5 ns. This pattern of spin polarization is in agreement with 
the experimental data of Fig. 4(c), Ref. ITU 

non-parabolicity of the energy spectrum. Remarkably, 
while this does not change the symmetry of the Hamil- 
tonian, the symmetry of responses is lower than in the 
special case (and in the naive picture). Our results give 
a microscopic explanation of experiments 11] and pro- 
vide a novel mechanism for generating spin polarization 
electrically via spin-orbit interaction. 

We consider a model of 2D electrons with charge e < 
and (pseudo-) spin i, obeying a Hamiltonian 

iI = efe-ib(k).o- + F(r), (1) 

where Ck is the dispersion law in the absence of spin-orbit 
couphng, V {y) is the potential due to impurites, V^(r), 
plus a small electric field E, cr are the Pauli spin ma- 
trices, and b (k) includes both intrinsic spin-orbit field 
bgo(k) and external field B. Wc consider in-plane mag- 
netic field, i.e., there is no orbital quantization, and dis- 
regard electron-electron interaction. In the following, we 
study spin polarization density s (r) = (cr) and spin 
currents (r) — {cr^, v})n2D- Here, ^20 is the elec- 
tron density, { , } is the anticommutator, and the velocity 
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•V = i [H, r] is spin-dependent. (We set ft = 1.) 

For a bulk 2D system with only intrinsic spin-orbit in- 
teraction, the kinetic equation has been derived 0, Q, 
ITsI l . Following Ref . Il5l we may write a spin-dependent 
Boltzmann equation for the distribution function, repre- 
sented as a 2 X 2 spin matrix / = /o(k) + i/c(k) 1 + 
f(k) • (T, with equilibrium distribution function /o, ex- 
cess particle density /c, k = {kx, ky) = (fccos(p, fcsin(^), 
and spin polarization density described by f . Magnetic 
field and spin-orbit coupling split the energy spectrum 
into two branches: for a given energy e, there are two 
Fermi surfaces. Thus, for elastic scattering, energy e is 
conserved but |k| is not, due to inter-branch scattering. 
In the following, we assume b <IC E-p. This motivates 
defining k^^ such that e^^ — e, and defining = e^. . For 
a fixed energy e, the velocity operator is |16| 



= k 



(1 + C) b • cr/2fc,] - (1/2) d{h ■ o-)/ak, (2) 



with unit vector k — k/fc and band nonparabolicity C, = 
{kf/v^){dve/dk^) — 1. Instead of using the distribution 
function / (k) as density in k-space, we consider it as 
a function of energy e and direction ip in k-space. In 
this representation, /c (k, e) and f (k, e) are transformed 
into distribution functions n{ip, e) and $ (</5, e), resp., 
which can be written as a matrix $ {ip, e); for a detailed 
derivation see Ref. The kinetic equation for E = Ex 
is [3 
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where fo is the Fermi distribution function, b = h{p) 
is evaluated for |k| = k^, and with charge distribution 
n = SPrVek^cosp and /3 = {eE/16n^Ve) {—dfo/de). In 
Eq. the first term is the partial time-derivative, the 
second term describes spin precession in the momentum 
dependent field b (p) , and the third term is the driving 
term, given in lowest order in E. 

The collision integral on the r.h.s. of Eq. (O can be 
found in Born approximation by Golden Rule jl5| , 



\ J coll 



dp' K{^) $((/3') - $((/3) 



+ / dp'<T-[M {p, n {ip') - M (^', ip) n (ip)] . (4) 
Jo 

Here, the first term describes spin-independent scatter- 
ing, with K{-d) ^ K {ip' -p) = W (q) kj2TTh^v,_ and 
q = 2fceSin(|i9|/2). The factor W{q) = {\Vi{q)\^) 
does not depend on the direction of the momentum 
transfer q since we assume that the system is macro- 
scopically isotropic. The second term in Eq. Q), de- 
scribed by Eq. (31) of Ref. [l^ includes two contribu- 



tions, arising from the spin-dependences of the density- 
of-states and of the momentum transfer for a fixed en- 
ergy e. These contributions are proportional to K {-d) and 
K {■&) = (dK/d-d) tan(i?/2), resp., and explicitly depend 
on p>, ip' through b {p) and b {ip'). 

We consider Rashba spin-orbit interaction, and choose 
the X axis along the field B, i.e.. 



b (k) = 2a z X k + Aj; X, 



.9*MbS, (5) 



with Zeeman splitting A^; , thus b (k) is in-plane and E 
and B are parallel, see Fig. 1(a). (For E — Ey there 
is yz mirror symmetry and Sz vanishes. Thus, the Sz 
term linear in E is determined only by the component 
Ex parallel to B.) 

Next, we write the kinetic equation Q in Fourier 
space by expanding the azimuthal dependence as / {ip) — 



Eoo 
m— - 



' fm ■ Combining the in-plane spin distribu- 



tion as $^((/?) + i^y{p}) = Eme'""^*m, and using the 
form of M given in Ref. we find |0| 

ilm = i^xKi + + (2 + 7(+)) Sm,2 

A^/3(l + 7o)(5|™|4 -T-ifc™^'™, (6) 

^* ) 



T-^kra^l,, (7) 



with inverse transport time t ^ — 2tt {Kq — Ki) and 
km = {Kq — Km)/{K(j — Ki). Also, 6q — 2iak, so the 



spin-orbit field can be written as bx + ib. 
Finally, the remaining parameters are 



y\B=o 



7(+) = 2(/c2 - 1), 7o = C + 2 



(8) 



In the limit of small- angle scattering, 70 = C + 3. We 
have assumed that B is time-independent and any time- 
dependence of E is slow compared to . 

Let us consider general properties of Eqs. ®-0- One 
can prove algebraically that 

'i>™-'i>^-„/ = = -{-!)"- (9) 

for both the stationary regime and for transients gener- 
ated by a time-dependent electric field. [Arbitrary initial 
conditions might deviate from Eq. @, but such devia- 
tions would decay to zero at least as fast as the spin re- 
laxation rate.] Also, these identities directly follow from 
the symmetry properties of the components of the pseu- 
dovector $ for the system with the axial symmetry Coov 
of the Rashba spin-orbit coupling in the fields E, B || x. 

In particular, the symmetry of Eq. © allows / 
in the stationary regime; therefore the spin polarization 
Sz is generally finite, despite the fact that the effective 
field b(k) has only in-plane components. Now we can 
evaluate Eq. for all m and Eq. ([7|l for m > 0, and 
eliminate complex conjugated quantities using Eq. (PJ. 
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Isotropic scattering, parabolic bands, stationary 
regime. First, we assume isotropic scattering and 
parabolic bands, thus Km — and 



that = ^T^y 



7(+),7o = 0, 



k — 1 — A 



(10) 



In this case, we solve the kinetic equations (jHl-IZll exactly 
by setting d^/dt = 0. The stationary solution is 

^-j^ = A^T/3<5|„|,i + 26„r/3 {6,n,o + S,n,2) , (H) 
= 0, (12) 

which can be checked by inspection. The total spin 
polarization density is s = s"'^ + s^, with equilibrium 
contribution 8°'^ || B and non-equilibrium contribution 
sj^ = 4tt J de^Q. Thus, the out-of-plane polarization van- 
ishes, Sz — 0, as one would expect from the above naive 
argument — even though the symmetry allows 7^ 0. 
Hence, vanishing s^ is a property of the specific model 
of Eq. (|10|l . On the other hand, even for this model, our 
solution is ^{f) — 4:t(3 [h{ip) — ■^A^x] cosip, i.e., in ad- 
dition to the total field b, there is a correction —^A^x, 
indicating that spin-orbit and external magnetic fields 
cannot be added. However, it does not contribute to the 
in-plane spin polarization, as it is averaged out when in- 
tegrating over (p. 

The polarization s°'', in the absence of spin-orbit cou- 
pling, arises from Pauli paramagnetism, s°'^ = ~ 
?i| = ^i^(e| — e|) = ^i^A^, with the density of states 

~ kp/nv-p — m*/7r, Fermi momentum kp, Fermi ve- 
locity vp, and effective mass m*. This spin polariza- 
tion does not depend on the electric field, thus <I>g = 0. 
On the other hand, the electric field causes drift, pro- 
ducing an average spin-orbit splitting, 6^^. = {by) ~ 
(1/'t-2d)// dedipn{e,ip)by = 2aeExT. By analogy to 
Pauli paramagnetism, one might guess that Sy = ^u{by). 
This expectation is indeed met, because Sy — aeE^Tv 
coincides with the value following from Eq. lllll . and it 
also agrees with known B — results 3, ^ Q, 0|- 
Hence, for the model of Eq. H1U|) . the in-plane polariza- 
tion can be described in terms of the average spin-orbit 
field. 

In the field A^,, the equilibrium spin polarization per 
electron is \uAx/n2P) — Ax/2Ep, so one expects that 
the drift caused by the charge current leads to a spin 
current = (j^/e) A^/^Ep. Our results agree with this 
expectation; evaluating the definition of by inserting v 
leads to = 2n J dev. Re [*i + + 2i3A^t] then 
Eqs. and (|ll|l are used. Note that the spin current 
is well defined for ba — because spin is conserved; 
also our calculations with finite 6^ result in the same j^. 
Other spin current components vanish, even for finite 
6a,; for Aj, = it is well-known that = [Hill [H 
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Anisotropic Scattering and Block equations. Now 
we consider anisotropic scattering and/or non-parabolic 
bands, and also include transients. We consider the 
"dirty limit," \ba\ ^ <C Ep with constant B such 



|A:r| < \ba\ < T- 



(13) 



where lo is the characteristic frequency of the field E, 
and Tz — ^Txy — i^ba^^r) are the Dyakonov-Perel spin 
relaxation times. In this regime, and ^'m decay ex- 
ponentially fast with increasing |to|, since 5^1^-^/ 5^ ~ 
r'^baAx ^ 1 for m > 2, and similarly for 5'^. This al- 
lows us to solve kinetic equations ©-0 order- by-order 
in the small parameter {t/tzY^'^. Considering the lowest 
non-vanishing order, it is sufficient to retain only equa- 
tions for |m| < 2. Eliminating the m = ±1, ±2 compo- 
nents yields the equations of motion for $ up to order 
(r/r,)V2 0. 

Finally, we evaluate the equations of motion for the 
total polarization s at low temperature T, taking all pa- 
rameters at the Fermi level. We obtain the Block equa- 
tion 



s = (b) X s - s + T, 



(14) 



where the spin relaxation tensor r ^ is diago- 
nal with components {T^y, T~y, t^^} and T = 
(ii^A^T-y\ \v{by)T-^, iz/Aj, (6y)7o). Note that our 
proof of Eq. H14(l is valid only in linear order in E [cf. 
Eq. (j3l], i.e., products {by) s^ were disregarded. 

To develop a physical picture for this central result, we 
note that Eq. (|14|l is a Bloch equation, where polarization 
s is generated with a rate F and then precesses in the to- 
tal field (b) = 5*/iBB + b(jr (Hanle effect). Most remark- 
ably, for anisotropic scattering and/or band nonparabol- 
icity, the combined effect of spin-orbit and external fields 
generates a spin polarization along the z axis with rate 
= \v g* ^lB.{^x{^s,o))zlo = \vaeExT A^^a, i.e., per- 
pendicular to botk magnetic and spin- orbit fields. This 
rate F^ arises as follows. Scattering of nonequilibrium 
carriers leads to an extra fc^-dependent x polarization 
due to the term proportional to 70 in Eq. On a 

timescale of r, this polarization then precesses around 
the y component of bso, as described by the first two 
terms of Eq. Q. 

Next we consider the dc case, s = 0. In the lowest 
order in E, the total spin polarization is s^ = ^vAx, 



Sy — —va cEt-t 
y 2 



1 



— va eEj.T — 

2 1-1- A^r T 



A^TxyTz 

70- 



70 



(15) 
(16) 



The first term of Eq. (|15|l arises from Eq. (|ll|l . while 
the second term and Sz are due to anisotropic scattering 
or nonparabolic bands. The dependence of Sz on A^; is 
in agreement with the data in Fig. Ic of Ref. |Oj where 
Ts = {TxyTzY'^'^ ~ 5 ns, suggesting that our microscopic 
model might explain the experimental observations. 

Spin currents. Evaluating Eq. © for m = 0, we 
find that $1;^ = —{Ax/ba)^Q. We then evaluate the 
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spin current at T = 0, jy = in f de v^lm^^^ = 

A-K / deve{l^xl\ba\)% = VY{^x/\ba\)sz, finding that 
is proportional to Sz- (This relationship also follows from 
the Heisenberg equation of dy.) Hence, the polarization 
Sz [Eq. ifTHjl ] leads to a transverse spin current j^; a fi- 
nite is in agreement with numerical results of Ref. |3 
For A.^ = 0, this relation is equivalent to the argument 
[21L l2a | based on equations of motion J^] , showing that 

Spin Dynamics. Even for isotropic scattering, a time- 
dependent electric field leads to an out-of-plane polar- 
ization Sz{uj) = ^uv/:^xh%{uj)/[l^l - + _ i[T-^ + 
T~^)ijj\] however, it has no static component Sz[uj = 0). 
Similar results were found in Ref. 25 for |6a| -C A^,, r^^. 

Spin dynamics is accessible in a pump-probe scheme 
[ni. Namely, spins can be pumped by applying a short 
electric pulse of duration tp <^ Tz, A~^. Then, accord- 
ing to Eq. the spin polarization immediately after 
the pulse is Sz{0) = tpTz oc Aa;7o, i.e., Sz(0) is an odd 
function of A^;. Solving the Bloch equation (|14|l . we get 



1 + 2/70 . ■ 

cos i It sm i It 



(17) 



with frequency fl — y^4A^^Tz)^^^/4Tz of Hanle oscil- 
lations (for consistency, we only consider terms linear in 
E). We plot Sz{t) in Fig. 1(b), taking the parameters of 
Ref.O and with a choice of 70 = 0.03, and find qualita- 
tive agreement with the experiment. The weak-field re- 
gion 4|Aa;|r2 < 1, where the oscillations are overdamped, 
is very narrow, \B\ < 0.25 niT. Note that the experimen- 
tal data shows that the sign of Sz depends on the sign 
of Ax, already on time scales much shorter than jAj,!^^. 
Therefore, the sign of Sz cannot be due to spin precession 
in the external magnetic field — implying that a polariza- 
tion generation mechanism like the one described above 
was experimentally observed in Ref. fill 

Strictly speaking, quantitative comparison with the 
data of Ref. ^3 cannot be performed because the films 
were of low mobility EpT 1, violating the assumptions 



of our Boltzmann description, and were in 3D regime 
(a coupling kyax — kxOy occurs here due to strain). 
Furthermore, in models with a more complicated spin- 
orbit interaction than the Rashba coupling, other sources 
of z— polarization might become important. However, 
Eq. (fT^ was satisfied, because h/xz '--^ 3 x 10^® eV; 
|A^j< 10-*^ eV; |6a| ~ lO'^ eV; and ft/r - 2 x 10^^ 

ev lillil. 

The effective field b(k) for a 2DEG with pure linear 
Dresselhaus coupling, on the (001) surface of a HI-V ma- 
terial, is obtained by replacing k on the right hand side of 
Eq. 10 by q = 7?.k, where 71 denotes reflection through 
the (110) crystal plane. Our result (|16|l for the polar- 
ization Sz can be applied to this case if we replace Ex 
by the component of the electric field along the direction 
B' = 7?,B. For general forms of the spin orbit coupling, 
we note that the C2V symmetry of the system ensures 
that if S = 0, there can be no term in Sz linear in E. 
However, there could be terms non-linear in E, if E is 
not parallel to a symmetry direction [110] or [llO], e.g. 
Sz oc E'^ — Ey where x refers to the [100] crystal axis, 
which would then give an all-electrical mechanism for 
generating out of plane spin polarization. 

In conclusion, we proposed a mechanism for generat- 
ing bulk spin populations polarized perpendicularly to 
magnetic and spin-orbit fields; for 2D systems this is an 
out-of-plane polarization. It relies on anisotropic impu- 
rity scattering and/or band nonparabolicity and provides 
a new method for electrical control of electron spins. Our 
model is derived for 2D systems, but the results should 
have a more general validity, and they agree with recent 
observations of combined effects of the external magnetic 
and spin-orbit fields in 3D samples. 
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I. LIST OF SYMBOLS 



Dispersion law in absence of spin-orbit interaction 
Vector of Pauli matrices, (cr^, ay, Oz) 
Total field in energy units, containing both spin- 
orbit and external magnetic fields 
Electric field, E = 
External magnetic field, B = _Bx 
Zeeman splitting, A^, = g* ^bB, with effective g- 
factor g* and Bohr magneton /iB 
charge of carrier, for electrons e < 
Rashba spin-orbit coupling constant, ba = 2iak 
Impurity potential 
Spin-dependent velocity, v = i [H, r] 
Anticommutator, {A, B} = AB + BA 
Spin polarization density, s (r) = 7120 (c) = s'^'' + 
s^, containing both equilirum and non-equilibrium 
contributions s^"* and s^, resp. 
Spin current, j'' (r) = n2n{\ {a^, v}). 
Wave vector. In two dimensions, k = {k^, ky) = 
fck = {k cos ip, k sin ip) 

Distribution function as 2 x 2 matrix, /(k) = 
/o(k) -f ^/c(k) 1 -I- f(k) ■ cr, as a function of wave 
vector k 

Equilibrium distribution function, spin-dependent 
due to magnetic field 
Non-equilibrium particle density 
Non-equilibrium spin polarization density 
Spin-independent wave number contribution for 
given energy e, i.e., eu, = e 

Spin-independent velocity contribution, Vt = e'f.^/h 
Band non-parabolicity, 1 -I- = (k, / v^){dvt / dk^) 
Density of states at Fermi level 
Non-equilibrium distribution function $ = 
$ (tp, e) as 2 X 2 matrix, as a function of direction 
of k and energy e, <1> = in 1 -|- $ ■ cr 
Excess particle density, as a function of ip and e, 
n = SPrv^k^ cos (fi 

Non-equilibrium spin polarization density, as a 
function of ip and e 

Angular dependence of spin-independent scat- 
tering, in Born approximation K (<^' — = 
W (q) k,/2-Kti'v, 
Scattering angle, = tp' — 
Momentum transfer, q = 2kt sin /2). 
Scattering contribution due to spin-dependence of 
momentum transfer, K (§) — {dK/d-ff) tan (i?/2) 
Fourier coefficients of $^(</?) 
Fourier coefficients of "1?^ (ip) -I- {ip) 
Fourier coefficients of K [■&) and K {■&), resp. 
Defined as fcm = {Ko-Km) / {Ko-Ki),i.e., km > 
for m > 

Transport lifetime, r^^ = 27r(A'o — Ki) 

Factor describing coupling to electric field, l3 = 

{eE/167T\,){-dfo/d€) 

Fermi distribution function 

Spin-dependent collision contribution due to b^, 
7(+) = 2(fe2 - 1) 



70 Spin-dependent collision contribution due to B, 
70 = C + 2 (A-i - Ko)/Ko - ifi 
c Coupling coefficient in kinetic equations, c = 
j6*A,r/3 

Tz Dyakonov-Perel spin relaxation time, Tz = 

Txy In-plane spin relaxation time, r^y = 2tz 
*r*7^ Spin relaxation tensor 

r Spin generation with rate F 



II. KINETIC EQUATION IN FOURIER SPACE 

In Fourier space, the kinetic equation Eq. (0) becomes, 
for m = 0, (note that fcg = 0) 



(SI) 



n=7^K^i-\Kn+c-^{'^^~n). (S2) 

2 z 2,1 



vifhile for rn ^ 0: 



^x/3o'5|m|,i. 



X (*„, - 



(S3) 

_ ^ 

" 2^ 



c 

2'5|m|,2, (S4) 



with — (3 + |/37(+) = ^2/? and with f3Q = f3 + (3jq. 

The contributions proportional to 7(_|_) and 70 arise 
from the second term in Eq. Q), where the kernel M is 
adopted from Ref. 0, 



keh{ip) 



b(v.) + b(^')^^^^ 



AhkfVe 



(S5) 



Concretely, we find 7(_|_) — 2'!it{K2 — Kq) and 70 = C + 
Attt{Ki~Kq)} By explicit evaluation of Km, we find the 
relations Kq = 2 X]„>o(~l)""^-^n and K„, = (-1)" A'q- 
\m\ Km ~ 2 J2o<n<\ni\ (-1)™^" uKn, m particular, Kq - 

ill = 2/?o+i^i and K2-Ka = 2{Ki-K2), see Sec UvTI 
This allows us to transform 7(.^) and 70 and we obtain 
Eq. ©. 



A. Solution for isotropic scattering 

We consider the stationary case d^m/dt ~ 0, isotropic 
scattering, and parabolic bands, thus 7(+) , 70 = and 



^ In terms of notation used in Ref. llSl we see that 7(-)-) = 7( 
for winding numer TV = 1 , corresponding to the field ba ; whereas 
70 = 7(+)/'?o |iv^]v=o fi'='<i ^- 
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Po = f3. We find the solution of the kinetic equation 



III. EFFECTIVE BLOCK EQUATIONS 



*m = A^T/3(5|„i|_i + 2baTf3 ((5™,o + Sm.,2) , (S6) 
'^L = 0. (S7) 



We now prove that Eqs. (|S6l)-(|S3 satisfy Eqs. (|gT|) - 
(IS4p . Eq. ljSl|l is trivially satisfied. Next, we use c = 
—ibaAxTf] and write 

*m = ^^|m|a + {^m,0 + <5m,2) , (S8) 



and insert into the l.h.s. of Eq. (|S2|I . 



i ic « , ic f 2ic 2ic\ 



.(i + i + l-l-l|c = 0. 



(SIO) 



For m ^ 0, we insert $f„ = in Eq. l(53(l . use /3(+) = fc2/3, 
and obtain 



In the following, we consider the dirty regime and solve 
the kinetic equations in orders of the small parameter 
(t/t^)^/^. The contributions to vj/^ and of order r" 

are denoted as and resp. It is convenient to 

choose units such that is of order unity. Let us now 
consider the regime, 



T-\ I A, I « «T-i, 



(S18) 



i.e., in Fourier space with respect to t, ^' 
—iuj'^m\uj) is of the same order in r as 

We take order 0{t-'^^^) of Eq. (|53|) for m = 2 

vI/(^/')-26„Tfc2-i/3(+), (S19) 

and order 0(r°) of Eq. HS4p for m = 1, 



(S20) 



r =26„fc2/3<^m,2 + A:r/3^|m|,l' (SH) 

which by comparison with Eq. ljS6p is also satisfied (be- 
cause fci = 1 by definition). Finally, evaluating Eq. (|S4() 
and dividing by c, we obtain 



0—^2 I *"|ni+l|,l + 



~2 (^*'^|ni-l|,l ~ ^ct^'^l™l,l 
~ f "^^"-S + T-'^-m.2 ) + Sun\.2 (S12) 

2i yAa; Aj; / 2 ' ' 

2'^|m+l|,l + 2'^|m-l|:l ^ '^"1.2 ~ <5m__2 + ^'^Im^Z 



•i/O 1 



i|,2 



(S13) 

= 0; 
(S14) 



therefore, our solution for <I> is valid. Note that the prop- 
erty kjn — ^ — Sm,o for isotropic scattering was not used 
explicitly in the above proof. 

Equations (|S6|) - (|S7| ) correspond to 

$2, {(fi) = 2t/3 {Ax - 2 |6q,| sinip) cos 93 = 4t/3 (cp) cos 93, 

(S15) 

<^y{<f) ^4:T\ba\f3 cos^ LP ^ At (3bl{ip) COS if, (S16) 



i^) = 0. 



(S17) 



where we have indicated above the equality sign that we 
used Eq. (|^ . Taking Eq. (EH for order 0{t^/'^) and 
using that = (^f)* yields 



(S22) 



(S23) 



iA,r^ 



z(l/2) ^-1 ,T,(l/2) 



*o ^ + 2r->a^/3, (S24) 



where we used T^y = ^r^ ^ = ^ \ba\^ t and — P- 

Next, we derive the equation of motion for <I>g. To this 
end, we take order 0(t") of Eq. HS3|I for to = 1, 



Then we take order ©(r^/^) of Eq. lf52)l . 



(S25) 



12)' 



(S26) 



25l I |2 *2(l/2) , * A I, \ 3 , 



A^ Im^-^ 



(1/2) 



l&cl' r$f - A, Imvl/(i/2) + c (2 + 70) , 



(S27) 

To) , 
(S28) 



As for the non-equilibrium charge distribution n, a fac- 
tor of cos(/3 is present. Remarkably, bs ((^) = b ((^) -f 
^ (7*/iBB diff'ers from the total field b [Eq. that en- 
ters in the Hamiltonian. 



^ The following derivation becomes simpler if we first insert the 
ansatz * = '^^^ + 5*. This replaces — > i^m\ — » S^^; 

— » 0; and /3o — > /37o in the following equations and c — > 
in Eq. J^6j. 
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where we used | A^jT (6* - ba) = ib'^A^r = c/p. 

Finally, using ib^ \b^fT^p = 2ib*^T~^Tl3 = 2cr-//A,, 
we obtain the Block equations for the nonequilibrium dis- 
tribution function 

^y(l/2) ^ ^.(1/2) _ 1 ^y(l/2) ^ I ,^ 

(S29) 
(S30) 

Therefore, the spin-orbit field generates y polarization 
with rate \ba\ t(3, whereas the combined effect of bso 
and B generates z polarization with rate c(2 + 7o), which 
is proportional to A^, . 



Sx = -T~y s.j; + ^i^A^ r-y\ (S38) 
Sy = A^s^- T-^ Sy + ^i^ {by) , (S39) 

Sz = {Sx - sf ) (fey) - A^ Sy - Sz + {by) A^ 7o. 

(S40) 

In linear order in E, {by) s^ vanishes and we can add the 
term s - { by) sf and {by) to the r.h.s. of Eqs. (|S38p 
and HS40|I , resp. Then, we can write the Bloch equations 
as 

s = (b) X s - s + r, (S41) 
with r = (i^/A, T-y\ \v {by) r-/, \u {by) A, 70). 



A. Bloch equations for total polarization s 

To obtain the spin polarzation s, we integrate the equa- 
tions for the total non-equilibrium contribution $ at tem- 
perature T = 0, see Sec. lIVBl Noting that the drift field 
is {by) = 2aeET, we find 



IV. PHYSICAL QUANTITIES EXPRESSED IN 
TERMS OF FOURIER COEFFICIENTS 



A. Transport lifetime 

The inverse transport lifetime is 



in f de2c—16Tr f deak^A^TB = — —aA^eET (S31) -1 f^^ T^ta\n q\ ^ i \ 

J J TTUF = / K {9) [1 - COS0) ^ 2Tr{KQ - Ki) . 

Jo 



= i^A^aeEr = {by) -vA,, = {by) (S32) 



Therefore, for the spin polarization l^pl j cj)^^ 



we obtain 



.bEB -IE 

'^x xy X 1 

^30l 



-1 „B 



+ 7^v{by)T^^, 



(S33) 
(S34) 



-v {by) A^7o, 



(S35) 



where we have used ib^A^ 



:7o 



^2ae£;rA,7o 



ji^ {by) Aa; 70- The stationary solution is 



^x'xyiz 



1 + Alr^yT, 
E 1 /L \ ^xTz 



■70 



(S36) 
(S37) 



We assume that A^, does not change over time, thus 
ds'^'^/dt — 0. The equilibrium polarization is Sy^ = s°'^ ~ 
and s%^ = \i'Ax, thus = -t~^s'^^ + \i>A^ T~y , which 



can then be added to the r.h.s. of Eq. (jS33p . This leads 
to the Bloch equations for the total spin polarization s = 



This motivated our definition 



Ko - Ki 

and one can see that km > for m > 1. 



(S42) 



(S43) 



B. Spin polarization density and spin currents 

We now evaluate the spin polarization density s (r) = 
I]a/3(^i (j")) ^pi^ Current (r) = 

Ea/3(V'i (r) H'^a/S' v}?A/3(r)}. The spin polarization 
density can readily be expressed in terms of the Fourier 
transformed distribution function, 



ded^ Tra^$(e, ^) = 4tt de $f„^o. (S44) 



For the spin current, one needs to evaluate 



= / / dedd Tr - {a„ v}<l(e, (S45) 



which is somewhat more complicated, as the velocity op- 
erator V depends on spin and on ip. It is obtained from 
the Heisenberg equation as 

■ rij- 1 FIT- « 1 ^'^fc 1 dh-(7 
V = i[H, x\ = - [H, dkj = — ^^71 — = vq + dv. 



dk 2 dk 



(S46) 
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For a fixed e, the value of the wave vector k depends on 
the spin, 



crb 

2v, 



(S47) 



[cf. Eq. (B35) in Ref.[T3. Thus, for vq = dek/dk = kvo 
and in lowest order in b. 



(S48) 

The velocity v at a fixed energy e can now be expanded 
in b, 



V = k 



(1 + C) 



b ■ (T 



2K 



+ (5v, 



(S49) 



yielding Eq. |j2Jl. [Equation is consistent with the 
gradient term ^[v- V/ + (V/) -v] in Eq. (B30) of Ref-E] 
Next, noting that 5 {c^, b • cr} = 6^^, we find 



v}ia,ki;, + k(l + C)^ 



2 dk' 



(S50) 



and we see that the second term in Eq. (jSSOp is not 
present in Ref. Il5l but would not lead to an extra con- 
tribution to j'' there. 

We insert Eq. HS50I) into Eq. (|S45p and use $ = in + 
cr • thus 



r = 



ded9 



2kw,$^ + k (1 + C)r^n- 



2h 



2 dk 



(S51) 
(S52) 



The components of the first term in Eq. (|S51|I are 



(S53) 
(S54) 
(S55) 
(S56) 



We evaluate Eqs. l)S54p and l)S54(l for concrete /.t, use 
'^m = + *'i'm' i&ke advantage of the fact that the 
spin current is a real quantity. This yields 



dede v^cosO^'^ie, 6) 
2^ J de«,($5' + $^J, 
2 JJ rfed6'wesin6l$^(e, 6*) 
2TTi I dew,($5' - $^1). 



■ly 



At: J deWeRe$i, 
-47r J de wjm$^. 



ii^'>'' ^2tt J de V, Re [$i + 



•'y 



2tt / de Im [$i - 



Jx 



2tt J de Ve Im [$i + , 
jWy ^2n [ de Re [$i - $_i] . 



(S57) 
(S58) 
(S59) 
(S60) 
(S61) 
(S62) 



Finally, we consider the remaining terms in Eq. I)S51|1 . Because (p = arctan(fcy/fca;), we have d(p/dkx 
and dip/dky = cos{if)/k, thus 



db. 



dk dh^ dip dh^ 



dkx dkx dk dk^ dip 



cos Lp 



db., 



dk dbfj_ dip dbfj 



sm p ■ 



db^ sin(/5 db^ 
'dk 
db 



fee dp} ' 

cospi dbfj. 



dky dky dk ' dky dp> " ^ dk k^ dp} 
We explicitly evaluate j*-^-*^ by inserting n = ( a/" ) cost/? into Eq. ljS51|l . 



,-(2)m 



•ly 



dede 



eET f^dfo 



47r2 V de 



db. 



db. 



cos^ p} (1 + O 6u — fee cos^ p} — + sin p cos (jS „ 

dk dp 



de ■ 



eET f ^dfo 
de 



8tt 



1 + C-k, 



d_ 

dk 



dede 



eET f-dfo 



47r2 V de 



{2b^,+bt^ + b''_,)-2{bt^ + b^_,) 
db, 



sin p cos p> {1 + b^j^ — kf sin pi cos pi -— p — cos^ p} 



dk 



dbj^ 

dp 



de ■ 



eET f -dfo 
de 



8tt 



i{-i + C-k, 



dk 



- shi{p>)/k 

(S63) 
(S64) 



(S65) 
(S66) 
(S67) 
(S68) 
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where we decomposed = e""'''6(^ and used cos^ ip — j {^2 + e _|_ g^"^) and sin 1^9 cosies = | (e 



-2iip _ p2iip\ 



Thus, for total spm currents j^, we find 



27T J de Ve l^i' + <P'l, +Pt (^1 + C - k,-^^ (26^ + 6^ + b'^a) - 2(6^ + ^^2) | 



For the model of Eq. 0, we then use 62 = ^'12 = ^ and 6q = Aj;(5^a;. 



(S69) 
(S70) 



For comparison, note that the charge current is Thus, we can use 

f JJ ded6 evTT^{e,e) = J J ded9 evn. (S71) 2^ (-1)" ni^„ = - 2 ^ {-i f nK^, (S82) 

n> |m| 0<n< |m| 



Evaluating the distribution /i, we get 



J. 



J J dede SefiTvlK cos^ e = j deSnePrv'^K (S72) 



and find 



K,^ = {-irKo-\m\K,^-2i-ir ^ {-i f nK^. 

T=0,C=0 2 El ^F^F e'^Er fcp e2i?Tn2D .o-70\ 0<n<|m| 

- ' ^^^-^2^-^^' (S83) 
i.e., we recover the Drude conductivity. 

C. Expressing Krn in terms of if„ 

Using that K [9] = J2n=-oo and that K^„, = 

Km, we get 

= tan ^ ^ in (e"^ - e"™^) X„ (S75) 

n>0 

= -2 ^ ni^„ sin (nO) tan -. (S76) 

n>0 

Because -RT (-0) ^ K {0), we get 

K,n = 7r ? ^ (^) 

Ztt J 2 

— —2 riKn I dO cos {mO) sin (n6') tan ■ 

^ — ' ZTT / 

„ ^ n *J 



n>0 



2 

(S78) 

For n > 0, we integrate 

Imn — — / M COS {mO) sin (n6') tan - (S79) 
271 j 2 

n < \m\ , 

4 n=|TO|, (S80) 

/ -|\l+n+m I I 

(-1) n> \7n\, 

which leads to 

km - \m\ Km + 2 (-1)™ (-1)" (^81) 



